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We systematically study the effect of high pressure on structure, electronic structure and ther- 
moelectric properties of 2H-M0S2, based on first-principles density functional calculations and the 
Boltzmann transport theory. Under pressure, cross-plane lattice size reduces much faster than in- 
plane one due to weak van der Waals interaction and it gets more difficult to shrink the size as ex- 
ternal pressure goes higher than 20 GPa, which agrees with experimental observation. A conversion 
from van der Waals bonding to covalent-like is found from charge density calculation. Concurrently, 
the dependence of band structure on pressure shows that a transition from semiconductor to metal 
occurs at 30 GPa. Band features close to Fermi level are discussed, such as narrow band from V to 
A and pressure-induced decrease of the band dispersion, which are advantageous for high values of 
thermopower. Our transport calculations also find that pressure-enhanced electrical conductivities, 
as well as high values of thermopower (up to a few hundred /iV/K), lead to significant values of 
thermoelectric figure of merit (above 0.10 for high pressure and even up to 0.65 at 30 GPa) over a 
wide temperature range. Our study supplies a new practical route to improve the thermoelectric 
performance of M0S2 and of other transition metal dichalcogenides by applying hydrostatic pressure. 

PACS numbers: 31. 15. A- 73.50.Lw, 72.80.Ga, 



I. INTRODUCTION 

M0S2, a typical member of transition- metal chalco- 
genide family, consists of alternating sandwiched sub- 
layers bonded by van-der- Waals (VDW) interaction. 
This weak inter-layer interaction, which has rendered it a 
good solid lubricanlfi, proves also essential in determining 
manifold novel properties, such as Raman frequency—'^, 
band-gap type (direct, indirect) and band-gap size^*^ 
changing with number of stacking layers. Besides the 
method of layer stacking, one also resorts to applying 
external pressure (stress) to further tune the inter-layer 
interaction^. To our best knowledge, pressure-induced 
effect on the VDW interaction has not been theoretically 
investigated. 

VDW interaction gives rise to some intriguing prop- 
erties on one hand, to a poor electrical transport 
properties^^^ on the other. M0S2 and other transi- 
tion metal dichalcogenide materials were reported to have 
very pronounced value of thermopower, but poor electri- 
cal conductivit}*2r"ii^, in result of a negligibly small value 
of dimensionless thermoelectric figure of merit (ZT ~ 
0.006 at the optimal doping and working temperature)^^. 
To improve upon its electrical conductivity for a better 
value of ZT, pressure/stress may be employed and ex- 
pected to suppress VDW bonding or to convert it to co- 
valent or metallic type. Encouragingly, similar studies 
have already been performed on the single layer. Stress- 
induced enhancement of electron mobility in M0S2 sin- 
gle layer has been demonstrated in experiments^ for a 
promising nano-electronic device. Both compressive — 
and tensile strain^^ in a single layer give rise to a re- 
duced band gap from first-principles calculations. A 



semiconductor-metal transition was obtained for a com- 
pressive strain of about 15% for single-layer 1^082^^. So 
far, there have been no reports on pressure effect on bulk 
electrical properties. Moreover, how pressure affects the 
electronic structure and if it can ultimately improve its 
thermoelectric properties have never been investigated. 

In this work, combining first-principles density func- 
tional calculations with semiclassical Boltzmann trans- 
port theory (BTT), we study the structure, electronic 
structure and thermoelectric transport properties of 
M0S2 under hydrostatic pressure. Our calculation shows 
a vanishing anisotropy of structural change at around 
20 GPa, agreeing with the experimental data^. We also 
find a concurrent conversion from VDW interaction to 
weak covalent-like type. A further compression induces 
a stronger covalent bonding and a semiconductor-metal 
transition at about 30 GPa. This transition improves 
greatly the inter-layer electrical conductivity but still 
keeps a large value of thermopower (~ 200/iy/i^), there- 
fore giving rise to much enhanced thermoelectric trans- 
port properties. Values of figure of merit (ZT) as high as 
0.15 for in-plane direction and as 0.65 along cross-plane 
direction are obtainable in a wide range of temperature 
(from 100 K to 700 K). Our study demonstrates that ap- 
plying high hydrostatic pressure can be an effective way 
to improve the thermoelectric performance of M0S2 and 
of other transition- metal dichalcogenides. 



II. METHODOLOGY 

M0S2 has Pea/mmc space group symmetry and con- 
sists of a hexagonal plane of Mo atoms sandwiched by 



2 





10 20 30 40 50 
Pressure (GPa) 



Mo-S B 
S-S • 

S1-S2 O 



[3-B-e □ □ B B- 



1.00 
0.95 
0.90 
0.85 
0.80 
0.75 
0.70 

I 

4.20 
3.70 
3.20 
2.70 
2.20 



10 20 30 40 50 10 20 30 40 50 
Pressure (GPa) Pressure (GPa) 



FIG. 1. (Color online) (a) Illustration of unit cell of M0S2 
and definition of the bond length (Imo-s, ds-s, dsi-S2 and 
angles and (j). The dependence on hydrostatic pressure of 
(b) lattice parameters a/ao and c/co, (c) bond angles, (d) 
bond lengths. A black dashed line in (b) is used to guide the 
eye. 



two hexagonal planes of S atoms. The unit cell contains 
two alternating layers with an AB stacking along c axis, 
as illustrated in Fig. [H 

Band structure of M0S2 is calculated by using the gen- 
eral potential linearized augmented plane- wave (LAP W) 
method as implemented in the WIEN2K package^^. The 
electronic exchange-correlation is described within the 
generalized gradient approximation (GGA) of Perdew- 
Burke-Ernzerhof(PBE) flavor which has been proved 
also efficient in dealing with the weak inter-layer inter- 
action within MoS2^. We use 5000 k points in the full 
Brillouin zone(BZ) to achieve a self-consistency accuracy 
better than 1 meV/atom. The Engel-Vosko GGA (EV- 
GGA) formalism^^ is applied to calculate the band gap 
more accurately. 

Transport properties was calculated based on the 
Boltzmann transport theory applied to the band struc- 
ture. The integration is done within the BOLTZTRAP 
transport code^^. A very dense mesh with up to 18000 
k points in the BZ is used. The electron scattering time 
is assumed to be independent of energy due to its good 
description of thermopower S(T) in a number of thermo- 
electric materials— 

We fully relaxed the unit cell volume, shape and inter- 
nal atomic coordinates under each hydrostatic pressure 
till the internal atomic force is less than 10~^ eV / A. We 
consider the external hydrostatic pressure no larger than 
50 GPa where no occurrence of structural phase tran- 
sition has been observed experimentally^ and structure 
remains to have Pea/mmc space group symmetry. 



III. RESULTS AND DISCUSSION 

We first studied pressure effect on structure of M0S2 
and understand the discontinuous structural variation 
observed in experiment^. Below 20 GPa, reduction of 
lattice parameters with pressure are anisotropic along 
in-plane and cross-plane directions, c/cq shrinks much 
faster than a/ao- This anisotropy vanishes above 20 GPa, 
which agrees with experiment by R. Aksoy et. al.^ and 
we analyzed by looking into internal structural parame- 
ters including bond angles and bond lengths, as shown in 
Fig. [Tl^c,d). In Fig. [TJc), opposite trends occur between 
the pressure dependencies of S-Mo-S bond angle and (j). 
Along with in-plane Mo-S and S-S bonds remaining rel- 
atively rigid in Fig.djd), pressure reduces the inter-layer 
distance via rotating Mo-S bond but very slightly. More 
importantly, d5'i_5'2 (the distance between adjacent sul- 
fur layers), presented in Fig. [TJd), shrinks abruptly with 
pressure especially below 20 GPa. This indicates a sub- 
stantial pressure induced change of inter-layer interac- 
tion. 

To see what happens to inter-layer interaction upon 
pressure, we studied total charge densities under hydro- 
static pressures ranging from to 50 GPa in Fig.[2l Along 
with d5'^_5'2 decreasing with pressure, VDW interaction 
seems to disapper, a weak covalent-like bond starts to 
form at 20 GPa and grows stronger. Compared to no 
charge distribution between two neighboring layers as 
marked in dashed white frame at GPa in Fig. Efa), 
electron charge becomes more and more converged be- 
tween the atoms Si and S2 as illustrated from Fig. [2fb) 
to (d). Obviously, VDW type bonding is converted to 
covalent bonding due to pressure. 

The change of inter-layer interaction will naturally lead 
to dependence of band structure on pressure from which 
we did find a transition from semiconductor to metal. 
Fig. [3] gives the calculated band structure evolving with 
pressure from to 50 GPa. Fermi energies of semicon- 




FIG. 2. (Color online) Total charge density p of M0S2 (a) at 
GPa, (b) at 20 GPa, (c) at 30 GPa and (d) at 50 GPa, seen 
from (110) cutting plane. A gradual change of bonding from 
(a) van-der-Waals type to (d) covalent-like type is seen. 
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FIG. 3. (a-e) EV-GGA electronic band structure of M0S2 
under different hydrostatic pressures. The size of band gap 
decreases gradually with external pressure till 30 GPa where 
a semiconductor to metal transition occurs, as illustrated in 
(f). Band gaps A and WrA are defined in (a). 



ductors are set at valence band maximum (VBM). The 
behavior of band gap is shown from (a) to (e) and sum- 
marized in (f). Indirect band gap A is initially between 
VBM at T point and conduction band minimum (CBM) 
at K point at zero pressure, but applying pressure will 
cause the CBM to shift to between K and F points. A 
drops linearly with pressure till a semiconductor-metal 
transition occurs at 30 GPa. This only takes a compres- 
sive strain of 7%, in contrast to 15% in single layer—. It 
seems much easier for the bulk material than the single 
layer to close the band gap A by applying pressure, due 
to the existing inter-layer interaction. 

It is worth noting that, an increase of band disper- 
sion with pressure is observed, except for two bands near 
Fermi level between F and A points. These two bands 
are actually zone-folded from one band as unit cell has 
one sub-layer. For convenience, both bands are treated 
as one with band width WpA, as defined in Fig. [3l^a). 
WpA decreases from 1 eV (at GPa) via 0.50 eV (at 
30 GPa) to 0.25 eV (at 50 GPa), an approximate lin- 
ear dependence of WpA on external pressure is seen from 
Fig.Hf). 

To understand the counterintuitive behavior of the 
band dispersion decreasing with pressure, we look into 
its electronic state at 30 GPa in Fig. SI In the energy 
range close to Fermi level, we are only concerned with 
Mo-d and S-p orbitals (Fig. Ill^a,c)), since they play pre- 
dominate roles in forming conduction channels for both 
in-plane and cross-plane directions. In Fig. |31^b,d), the 
two cross-plane bands with the anomalous width WrA 
are due to hybridization of Mo-d^^ and S-p^ orbitals. 
The valence band at F point has a F4 symmetry -5^*2^, 
and the band below has a dominating contribution of 
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FIG. 4. (Color online) Electronic band character at hydro- 
static pressure of 30 GPa. The electronic states close to Ef is 
mainly from (a) Mo-d and (c) S-p orbitals. (b) S-p^ and (d) 
Mo-d22 play major roles for cross-plane conduction channel. 
Possible combinations of S-p^ and Mo-d^2 orbitals in M0S2 
unit cell is illustrated in (e) with symmetry labels used for 
the irreducible representations of the PGs/mmc space group 
(Ref. 22). 



S-p^ orbitals with F^ symmetry2^i2i at F point. The 
symmetries of possible combinations of Mo-d22 and S- 
Pz orbitals are illustrated in Fig. Hl^e). We found the 
valence band have S-p^ anti-bonding characteristic. In 
this sense, the shrinkage of inter-layer distance d5'i_5'2 
with pressure suggests a stronger anti-bonding state and 
therefore a narrower band dispersion near from F to 
A. It also indicates that the band with F^ symmetry at F 
point contributes to bonding. As pressure increases from 
30 GPa, the Fermi level shifts down from VBM by un- 
populating the antibonding state between F and A and 
is expected to strengthen the bonding state between sub- 
layers, which agrees with the charge density dependence 
on pressure as depicted in Fig. [2l 

From the band structures shown above, we can also 
extract some useful information on the electrical conduc- 
tivity. In Fig. [3fc-e), we found Fermi velocity (vi?) and 
effective mass (m*) respectively increases and decreases 
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FIG. 5. (Color online) Dependence of thermopwer on the 
temperature under different hydrostatic pressures along (a) 
in-plane direction and (b) cross-plane direction, respectively. 



with pressure. Calculated electronic densities of states 
(DOSs) (not shown here) close to Fermi level at 30, 40 
and 50 GPa are 0.07, 0.30 and 0.60 electrons/eV/u.c, 
respectively. Under thermal excitation within a range of 
interest (e.g., 700 K, equivalent to 0.06 eV of energy), 
DOSs translate to carrier concentrations n of approxi- 
mately 5x10^^, 2x10^^ and 5x10^^ m"^ for 30, 40 and 
50 GPa, respectively. The boost of carrier concentration 
(n) by pressure, together with increased v^ and reduced 
m* , suggest pressure-enhanced electrical conductivity. 

Also information on thermopower can be extracted 
from the band structures. According to the Mott 
relation^^, thermopower depends crucially on the deriva- 
tive of logarithmic electrical conductivity over energy 
close to Fermi level^^^^'^, namely, S ^^^\e=Ef- Based 
upon this, at the transitional pressure, a high value 
of thermopower is anticipated for both in-plane and 
cross-plane directions. Moreover, the anomalous bands 
upon pressure get more localized due to inter-layer anti- 
bonding states existing and its band edges have a big 
chance to touch Fermi level. It is therefore more ad- 
vantageous for having high values of thermopower along 
cross-plane direction. Meanwhile, the carrier concentra- 
tion (n ~ 10^^ m~^) obtained above is lower by orders 
of magnitude than that in good metals (n ~ 10^^ m~^). 
The lower level of carrier concentration, the higher value 
of S^^. Much higher values of thermopower will be ex- 
pected in our metallic M0S2 with pressure above 30 GPa 
than that in good metals, which are actually confirmed 
by our following calculations. 

We discussed qualitatively how pressure affects ther- 
mopower and electrical conductivity from band struc- 
ture. Now we turn to semi-classical Boltzmann transport 
theory to quantify thermoelectric parameters of M0S2. 

Firstly, we calculated thermopower as function of pres- 
sure and temperature. In Fig. O high values of ther- 
mopower were indeed obtained from 30 to 50 GPa. Es- 
pecially, over 100 /iV/K in the plane and more than 200 
/iV/K along c-axis are achievable over a wide tempera- 
ture range in the region close to 30 GPa. Also from Fig.[5l 
a strong anisotropy of thermopower between two per- 
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FIG. 6. (Color online) (a) The ratio of electrical conductivity 
cr over electronic relaxation time r and (b) -p^ as function 
of pressure and temperature. (50 : xx) used in (a,b) means 
(pressure : direction) . The ratio of Lorenz number Lo over ^ 
along (c) in-plane and (d) cross-plane direction under various 
pressures above 30 GPa is given to test if the Wiedemann- 
Franz law fits here. 



pendicular directions appears when pressure increases. 
The in-plane thermopower drops substantially while the 
cross-plane thermopower goes down fairly modestly. 

In Fig. ini^a), we then show calculated a/r as function 
of pressure and temperature, a/r increases with external 
pressure, so does cr if r is approximately pressure inde- 
pendent. To see how temperature dependence of a but 
without available experimental r, we assume electrons 
are mainly scattered by phonons at high temperature and 
electronic relaxation time r ~ T~^ is used^^^^^. Hence 
cr/(rT) as function of temperature can give information 
on temperature dependence of a. We show (j/(rT) ver- 
sus temperature in Fig. [Sfb). Typical metallic behavior 
is seen by an inverse power law relation between (j/{tT) 
and temperature, except for at 30 GPa where a slight 
deviation from the inverse power law is observed and 
suggests its conductive properties different from typical 
metals. 

Electronic thermal conductivity He/r was also calcu- 
lated but not given here. Alternatively, we put and 
cr in the background of the Wiedemann-Franz law and 
draw a ratio of Lorenz number (Lq = 2.8 x 10~^ Wl^/K^) 
over L (= ^) as function of temperature and pressure 
in Fig. [6l^c,d). The ratio shows a better obedience of the 
Wiedemann-Franz law as pressure increases from 30 to 
50 GPa, Lo/L approaches unity as pressure gets to 50 
GPa. In contrast to anisotropic thermopower S for two 
perpendicular directions, isotropic behaviors are found of 
dependence of L on temperature and pressure between 
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FIG. 7. (Color online) Temperature and pressure dependence 
of figure of merit ZT of conductive M0S2 along (a) in-plane 
and (b) cross-plane directions. The cross-plane is more prefer- 
able for practical thermoelectric application. 



in-plane and cross-plane directions. 



As is known, dimensionless thermoelectric figure of 
merit ZT is defined as ^ (or 7^), with tii as the 
thermal conductivity of lattice part and ^ = As- 
suming Ki is insignificant to Ke for metallic state^^ above 
30 GPa, we show calculated ZT value in Fig. [71 Value 
of ZT can reach 0.15 along in-plane direction and 0.65 
along c axis over a wide range of temperature at 30 GPa. 
As pressure increases, the values of ZT decrease. The in- 
plane ZT drops very abruptly, while the cross-plane ZT 
still has value higher than 0.10 over a wide range of pres- 
sure, suggesting a proper range of pressure is essential to 
optimize the thermoelectric properties. 



IV. CONCLUSION 

In summary, we systematically study the structure, 
electronic structure and thermoelectric properties of 2H- 
M0S2 under hydrostatic pressure, based on the Boltz- 
mann transport theory and first-principles density func- 
tional calculations. We find it becomes more difficult to 
reduce size along the cross-plane direction under pres- 
sure higher than 20 GPa, because of a crossover from 
weak inter-layer VDW interaction to covalent-like bond- 
ing. A quantitative agreement between our calculations 
and experimental data is reached. A semiconductor- 
metal transition is found to occur at 30 GPa and a much 
smaller compressive strain in bulk material than in sin- 
gle layer is needed for closing the band gap. Along with 
the pressured-enhanced electrical conductivity, we also 
obtained very large values of the thermopower (up to a 
few hundred /iV/K from 100 to 700 K) and ascribed it 
to some anomalous band features. Above all, good val- 
ues of figure of merit are found. Besides, an anisotropy 
of ZT value is found to exist between two perpendicular 
directions. The cross-plane direction is more preferable. 
Its ZT value under pressure can be obtained as big as 
0.10 and even up to 0.65 at 30 GPa over a wide temper- 
ature range. Our study supplies a new practical route 
to improve the thermoelectric performance of M0S2 and 
of other transition metal dichalcogenides by applying hy- 
drostatic pressure. 
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significant here than electronic Ke, due to following reasons. 
(1) At low temperature, electronic thermal conductivity Ke 
dominates in the metals. (2) At high temperature, Ke is 
insensitive to temperature under the circumstance of r ^ 
T~^, while Ki decreases with temperature, as experiments 
suggested^°. It makes even more significant than at 
high temperture. (3) It would reduce our estimation of ZT 
by only one half if ki were comparable to Ke. 



